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ABSTRACT 

We determine the space density of UV-luminous star-burst galaxies at z ~ 6 using deep HST 
ACS SDSS-i' (F775W) and SDSS-z' (F850LP) and VLT ISAAC J and band imaging of 
the Chandra Deep Field South. We find 8 galaxies and one star with [i' — z') > 1.5 to a depth 
of 

^AB — 25.6 (an 8 a detection in each of the 3 available ACS epochs). This corresponds to 



an unobscured star formation rate of 
Lyman break population at z = 3 — 4 (Ha 



15/ifo"M0yr- 



0.7, n 



M 



at z = 5.9, equivalent to L* for the 
= 0.3). We are sensitive to star form- 



ing galaxies at 5.6 < z < 7.0 with an effective comoving volume of « 1.8 x 10^ h^Q Mpc^ 
after accounting for incompleteness at the higher redshifts due to luminosity bias. This vol- 
ume should encompass the primeval sub-galactic scale fragments of the progenitors of about 
a thousand L* galaxies at the current epoch. We determine a volume-averaged global star for- 
mation rate of (6.7 ± 2.7) x 10^** hjo Mq yr~^ Mpc^^ at z ~ 6 from rest-frame UV selected 
star-bursts at the bright end of the luminosity function: this is a lower limit because of dust 
obscuration and galaxies below our sensitivity limit. This measurement shows that at z 6 
the star formation density at the bright end is a factor of ^ 6 times less than that determined 
by Steidel et al. (1999) for a comparable sample of UV selected galaxies at z = 3 — 4, and 
so extends our knowledge of the star formation history of the Universe to earlier times than 
previous work and into the epoch where reionization may have occurred. 

Key words: galaxies: formation - galaxies: evolution - galaxies: starburst - galaxies: high 
redshift - ultraviolet: galaxies - surveys 



1 INTRODUCTION 

How, when, and over what time-scale galaxies formed are questions 
at the forefront of both observational and theoretical cosmology. In 
recent years there has been tremendous progress in the study of 
galaxy formation from an observational perspective. Spectroscopic 
samples of star-forming galaxies have been detected at progres- 
sively higher redshifts: z ~ 1 (Cowie et al. 1996); 2 ~ 3 (Steidel 
et al. 1996); 2 ~ 4 (Hu & McMahon 1996; Steidel et al. 1999); 
beyond 2 = 5 (e.g., Spinrad et al. 1998; Weymann et al. 1998; Hu, 
McMahon & Cowie 1999; Dawson et al. 2001; Lehnert & Bremer 
2003) and now up to 2 = 6.6 (Hu et al. 2002; Kodaira et al. 2003). 

Complementing this ground-based spectroscopic work there 
have been a number of programmes based on the Lyman break 
"photometric redshift" technique (e.g., Steidel, Pettini & Hamil- 
ton 1995). Hubble Space Telescope (HST) imaging with WPC2 
and NICMOS has extended such studies to fainter magnitudes 
than accessible through ground-based spectroscopy (e.g., Madau 
et al. 1996; Dickinson et al. 2000; Thompson, Weymann & Storrie- 
Lombardi 2001). From the results of this work, great progress has 



been made in reconstructing "the star formation history of the Uni- 
verse" - the evolution of the comoving volume-averaged star for- 
mation rate (as shown on the 'Madau-Lilly' diagram: Lilly et al. 
1996; Madau et al. 1996) over the redshift range — 4. 

Here we describe a survey that extends such work to 2 ~ 6. A 
survey may be considered useful for selecting candidate star form- 
ing galaxies if it satisfies the following primary criteria: (i) it can 
distinguish 2 ~ 6 galaxies from lower redshift objects; (ii) it is 
deep enough to reach star formation rates characteristic of this pop- 
ulation; (iii) the volume is sufficient to encompass a statistically 
meaningful sample. 

In this paper we show that public data from the Advanced 
Camera for Surveys (ACS, Ford et al. 2002) on HST, released 
as part of the public Great Observatories Origins Deep Survey 
(GOODS - Dickinson & Giavalisco 2002^) programme, is of suf- 
ficient depth and volume to detect very high redshift galaxies, and 
that colour selection should reject less distant objects. We use this 
data set to examine the space density of UV-luminous starburst 
galaxies at 2 ~ 6. Observational evidence of the Gunn-Peterson 
effect in 2 ~ 6 QSOs (Becker et al. 2001) suggests that this is the 
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epoch of reionization, where starlight may for the first time domi- 
nate the Universe - although recent results on the polarization of the 
cosmic microwave background from the WMAP satellite indicate 
that the onset of ionization may have been at even higher redshift 
(Kogut et al. 2003). Using deep imaging by ACS through its SDSS- 
i' and -z' filters, we identify a sample of objects which are likely 
to be Lyman break star forming galaxies at 5.6 < z < 7.0. From 
the UV continuum longward of Arcst = 1216 A detected in the z' - 
band, we place a lower limit on the global star formation rate at 
2 ~ 6 and hence extend the Madau-Lilly diagram to early epochs. 

Throughout we adopt the now widely accepted 'concordance' 
cosmology: a A-dominated, flat Universe with SIa = 0.7, Q.m = 
0.3 and Hq = 70 /irokms"^ Mpc~\ The HST/ACS magnitudes 
in this paper are quoted in the AB system (Oke & Gunn 1983). 



2 OBSERVATIONS 

2.1 The ACS data 

We have analysed the first three epochs of GOODS data from the 
Chandra Deep Field South (CDFS) taken with HST/ACS^. Each 
epoch of data comprises a half-orbit i'-band image taken through 
the F775W filter (1040 sec, split into 2 for cosmic ray rejection) 
and a full-orbit z'-band image observed with the F850LP filter 
(2010 sec, split into 4). We use the vO.5 release of the reduced data. 

The released data set has been pipeline-processed for dark 
and bias subtraction and flat-fielding, and known hot pixels, bad 
columns and other cosmetic defects flagged in a data quality file. 
The post-pipeline images are sky-subtracted,"drizzled"(Fruchter & 
Hook 2002) and corrected for geometric distortion. Cosmic rays 
were also identified and masked. 

We use the zero-points for AB magnitudes deter- 
mined by the GOODS team: mag^g — zeropoint — 
2.5 logjo(Count rate/s~^) where the zeropoints for i'-band 
(F775W) and z'-band (F850LP) are 25.66 and 24.92 respectively. 
We have corrected for the small amount of foreground Galactic 
extinction toward the CDFS using the COBE/DIRBE & IRAS/ISSA 
dust maps of Schlegel, Finkbeiner & Davis (1998). The optical 
reddening is E{B — V) — 0.008, equivalent to extinctions of 
Af775 = 0.01 7 & Afsfmlp = 0.012. 

2.2 Construction of the catalogues 

Candidate selection for all objects in the field was performed us- 
ing version 2.2.2 of the SExtractor photometry package (Bertin & 
Amouts 1996). As we are searching specifically for objects which 
are only securely detected in z', with little or no flux in the i'- 
band, fixed circular apertures I'.'O in diameter were trained in the 
z'-image and the identified apertures were used to measure the 
flux at the same spatial location in the i'-band image, running 
SExtractor in two-image mode. For object identification, we de- 
manded at least 5 contiguous pixels above a threshold of 2a per 
pixel (0.01 counts pixel"^ s"'^) on the drizzled data (with a pixel 
scale of O'.'OS pixel"^). From the output of SExtractor, we create a 
sub-catalogue of all objects brighter than z^b < 25.6 mag, which 
corresponds to an 8 a detection in each epoch. 

In the compilation of the catalogues a number of 'figure-eight' 
shaped optical reflections and diffraction spikes near bright stars 
were masked, as was the gap between the two ACS CCDs and 



its vicinity in each tile. This means that the survey area is non- 
contiguous. 

Excluding the masked regions, the total survey area in the 
15 tiles of CDFS is 146arcmin^. In this region, we detect 10728 
sources in the z'-band brighter than z^b = 25.6 (our 8 a cut for 
the F850LP image from each individual epoch). The IRAF 'Art- 
data' package was used to create artificial objects of known magni- 
tude and we confirmed that we recover >95 per cent of such objects 
at the magnitude limit of our catalogue. 



3 CANDIDATE Z > 5.6 GALAXIES 

3.1 Redshift Discrimination 

In order to select z > 6 galaxies, we use the Lyman break technique 
pioneered at z ~ 3 from the ground by Steidel and co-workers and 
using HSThy Madau et al. (1996). The usual Lyman break tech- 
nique at z ~ 3 — 4 involves three filters: one below the Lyman 
limit (A.ost = 912 A), one in the Lyman forest region and a third 
longward of the Lyman-a line (Arcst ~ 1216 A). At higher red- 
shifts, we use just two filters, since at z ~ 6 the integrated optical 
depth of the Lyman-a forest is ^ 1. Hence the continuum break 
at the wavelength of Lyman-a is large, as is illustrated in FigureQ 
and the shortest-wavelength filter below the Lyman limit becomes 
redundant. Thus provided one works at a sufficiently-high signal- 
to-noise ratio, i'-band drop-outs can be safely identified through 
detection in a single redder band (i.e., SDSS-z'). This has been el- 
egantly shown by Fan and the SDSS collaboration in the detection 
of z ~ 6 quasars using the i' - and z'-bands alone (Fan et al. 2001). 
The sharp sides of the SDSS filters help in clean selection using the 
photometric redshift technique. 

In Figures |2|&|3| we illustrate how the use of a colour cut of 
(i' — z')ab > 1.5 can be used to select objects with z > 5.6, and 
we have 9 at z^b < 25.6. Objects satisfying our selection criteria 
in any of the three epochs were gathered into a preliminary cata- 
logue. The flux of these objects in the three epochs was averaged 
for each waveband, and those objects still satisfying the criteria 
Zab < 25.6 and (i' — z')ab > 1.5 in the combined data set were 
selected. A number of objects appearing in catalogs generated by 
the SExtractor software were identified by eye as different regions 
of the same extended source (visible in the ACS F606W t;-band 
or F775W i'-band images) and eliminated from the final selection. 
The 9 remaining objects satisfying our selection criteria are detailed 
in Table 1 . 

3.2 Interlopers at Different Redshifts 

Using our (i' — z') colour-cut criterion alone may also select ellip- 
tical galaxies at z ~ 2 (as shown by Figure|2}, the Extremely Red 
Object (ERO) population (e.g., Cimatti et al. 2002). Another poten- 
tial source of contamination are Galactic L- or T-dwarf low-mass 
stars (e.g., Hawley et al. 2002). We can guard against these lower- 
redshift interlopers by considering the near-infrared colours. Deep 
J and Ks band observations of the ACS survey field have been ob- 
tained with ISAAC on VLT, and reduced data from this ESO Imag- 
ing Survey has also been publically released as part of the GOODS 
programme^. As Figure |4| illustrates near-infrared colours effec- 
tively discriminate the ERO population from high redshift objects. 
Of the five candidates satisfying our selection criteria for which 
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Figure 1. The ACS-i' and -z' bandpasses overplotted on the spectium of a 
generic 2 = 6 galaxy (soUd hne), illusti'ating the Lyman break technique. 



such data is available, one (object 2) is detected in F606W t;-band 
(vab = 27.0 ±0, .2) and has colours consistent with being an ERO 
contaminant. All other candidates have vab > 28.0 (our 3 a detec- 
tion limit), although object 9 also has red near-IR colours. Object 
5 has a colours which suggests that it is likely to be a low mass 
star: (^ab ~ Jvcga)=3.12 & (J - A's)vcga = 0.2 (an L- or T- 
dwarf)- Comparison with the sample in Hawley et al. (2002) indi- 
cates a spectral type around T3. However, of the remaining objects 
for which no reduced infrared data is yet available, one is now con- 
firmed as a high redshift object (see 'Note Added in Proof) and the 
others are all fully resolved as shown in Figure|5| The contamina- 
tion rate would thus appear to be about 25 per cent. 



3.3 Limiting Star Formation Rates 

We are sensitive to objects above z ~ 5.6 and in principle could 
identify objects out to 2 ~ 7.0 where the Lyman-Q line leaves the 
z' filter. The rapid falloff in sensitivity of the ACS detectors above 
9000 A, however, combined with the transition of the Lyman-a 
break through the z' filter and the resulting incomplete coverage of 
that filter, means our sensitivity to star formation drops rapidly past 
z « 6.0, and renders us unlikely to detect anything past z « 6.5 
(Figure|6). The effect of luminosity- weighting of the redshift range 
also affects our effective survey volume, and we quantify this in 
SectionIO 

The observed luminosity function of Lyman break galaxies 
(LBGs) around Arcst = 1500 A is = 24.48 at (z) = 3.04 and 
m*j = 24.48 at (z) = 4.13, from Steidel et al. (1999), with a faint 
end slope a = —1.6 and normalization 0* ~ 0.005 /170 Mpc~^. 
If there is no evolution in the luminosity function from z = 3 (as is 
found to be the case at 2 = 4) this would predict an apparent mag- 
nitude of = 25.6 for an Llbg at 2 ~ 6. Thus our complete 
catalogue of z'^-g < 25.6 would include galaxies down to Llbg at 
2 ~ 6 if there is no evolution in the LBG luminosity function. 

The relation between the flux density in the rest-UV around 
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Figure 2. Model colour-redshift tracks for galaxies with non-evolving stel- 
lar populations (from Coleman, Wu & Weedman 1980 template spectra). 
The 'hump' m(i' — z') colour seen at ^ 1 — 2 is due to the 4000 A break 
redshifting beyond the i'-filter. 
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Figure 3. A colour-magnitude diagram for our data with the 2 '-band limit 
of < 25.6 and the colour-cut of (i' — z')ab = 1-5 shown (dashed 
lines). Candidate objects are marked with circles, or 3 cr lower-limits on 
the (j' — 2') colour where objects are undetected in i' . Note: A catalogue 
based on a simple colour cut is contaminated by stars, EROs and wrongly 
identified extended objects as described in the text. 
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Figure 4. Model (2' — J) - (J — K) tracks for galaxies with non-evolving 
stellar populations. Redshift intervals are marked, and we plot the colours 
of our (i' — z')ab > 1-5 objects for which near-infrared data is available. 
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Figure 5. A magnitude-size diagram for our data illustrating the clearly 
defined stellar locus at Rf^ < 0.1". Candidate objects are marked with 
circles — 6 are clearly resolved, and 3 are compact (of which one is likely 
to be a T-dwarf star, and another a putative AGN). 



Figure 6. The limiting star formation rates reached as a function of red- 
shift for our catalogue with < 25.6 mag. The star formation rates 
are inferred from the rest-frame UV continuum around 1500 A (Madau, 
Pozzetti & Dickinson 1998). We account for the filter transmission, and the 
blanketting effects of the intervening Lyman-a forest. We introduce small 
fc-corrections from the observed rest-wavelengths redshifted into the z'- 
band to 1500 A: the solid line assumes a spectral slope (3 = —2.0 (where 
f\ oc A*^ ) which is appropriate for an unobscured starburst, and the dotted 
line has (3 = —1.1 (the slope for mean reddening of the 2 ?a 3 Lyman 
break galaxies). The value of SFR^y = 15 h^^ Mq yr~^ for 2 = 3 — 4 
from Steidel et al. (1999) is the horizontal line. Our colour selection should 
remove most 2 < 5.6 galaxies (solid vertical line), and our effective survey 
volume extends out to 2 = 6.1 (dashed line). 



ilar mass range. From Steidel et al. (1999), the characteristic star 
formation rate at the knee of the luminosity function, LlbG' is 
SFRuv = 15 Mq yr~^ for z — 3—4. Figure|6|shows our lim- 
iting star formation rate as a function of redshift, inferred from the 
rest-frame UV flux in the z'-filter. We account for the filter trans- 
mission, and the blanketting effects of the intervening Lyman-a 
forest. We introduce small fc-corrections to Arcst = 1500 A from 
the observed rest-wavelengths longward of Lyman-a redshifted 
into the z'-band: we consider both a spectral slope (3 — —2.0 
(where fx oc A'') which is appropriate for an unobscured starburst 
(flat in /y), and also a redder slope of /3 = —1.1 (approriate for 
mean reddening of the 2 « 3 Lyman break galaxies, Meurer et al. 
1997). At limiting magnitude of z^b — 25.6, we can detect un- 
obscured star formation rates as low as 15 [16.5] h^Q Mq yr~^ at 
5.6 < 2 < 5.8 and 21 [25] /i^ Mq yr"^ at z < 6.1 for spectral 
slope P = —2.0 [—1.1]. Our (i' — z')ab > 1-5 colour cut should 
remove most galaxies at 2: < 5.6. 



1500 A and the star formation rate (SFR in Mq yr ^) is given 
by Luv = 8 X lO^^'SFRergss"^ Hz"^ from Madau, Pozzetti & 
Dickinson (1998) for a Salpeter (1955) stellar initial mass func- 
tion (IMF) with 0.1 < M* < 125 Af©. This is compara- 
ble to the relation derived from the models of Leitherer & Heck- 
man (1995). However, if the Scalo (1986) IMF is used, the in- 
ferred star formation rates are a factor of ~ 2.5 higher for a sim- 



3.4 Candidate Objects 

Using the selection criteria described above we have identified 7 
candidate galaxies photometrically selected to lie at 2 > 5.6. These 
(along with the probable star and ERO included for comparison) are 
described in Table 1 and cut-out images of each object in the 2', i'. 
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Figure 7. Candidate z > 5.6 galaxies selected to have (i' — z')ab > 1-5. 
The images shown are 3 arcseconds to a side and represent data from all 
three epochs of observation. Note that object 6 is not observed in epochs 1 
and 3 due to the nature of the tiling scheme used. 



«(F606W) and (where available) J and Kg bands are presented in 
FigureQ 

As Figure |5|illustrates most of our (i' — z')ab > 1-5 objects 
are spatially resolved in the HSTI ACS images (6 of the 9). Of the 
three compact objects, object 5 has the colours of an L- or T-dwarf 
star. We note that object 4 is almost 1 magnitude brighter than the 
next brightest object and lies close to the FWHM stellar locus, and 
is thus likely to be an AGN. 



4 GLOBAL STAR FORMATION AT Z ~ 6 
4.1 The Survey Volume 

Galaxies in the range 5.6 < z < 7.0 would be selected by our 
colour-cut provided they are sufficiently luminous. However, in 
practice we are sensitive over a smaller range of redshifts, because 
galaxies become progressively fainter at higher redshift, and so an 
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Figure 8. The 'Madau-Lilly' plot illustrating the evolution of the comoving 
volume-averaged star formation rate, reproduced from Steidel et al. (1999) 
and recalculated for our cosmology and higher luminosity limit of 1 L* (in- 
stead of 0.1 L*). The luminosity functions for galaxies derived by Steidel 
et al. were used to recalculate these values, assuming a slope of o = —1.6 
for z > 2 and a = —1.3 for z < 2. Data from the CFRS survey of Lilly 
et al. (1996) are shown as open circles; data from Connolly et al. (1997) are 
squares; the HDF-North results from Madau et al. (1996) are triangles; and 
the Lyman break galaxy work of Steidel et al. (1999) is plotted as crosses. 
The result from this work is shown as a filled circle. Note that none of the 
data-points on this figure are corrected for extinction. 



ever smaller fraction of them are able to satisfy the magnitude limit 
of z'j^-g — 25.6. There are three main reasons for this: 

(i) The selected galaxies lie primarily on the steep exponen- 
tial cut-off of the luminosity function (brighter than L* for a no- 
evolution model), so small changes in the absolute magnitude limit 
with redshift will greatly affect the number of galaxies brighter than 
the limit. 

(ii) Another effect is the fc-correction (Section [3.3> : as redshift 
increases, the z'-band samples light in the rest-frame of the galaxies 
at wavelengths that are increasingly far to the blue of 1500 A, where 
the LBGs' luminosity function was calculated. This fc-correction 
makes higher redshift galaxies fainter (as LBGs are red, with (3 > 
— 2.0), and therefore less likely to satisfy the magnitude limit. 

(iii) At redshifts z > 6, Lyman-a absorption from the forest en- 
ters the z'-band and makes galaxies fainter still, as there is incom- 
plete coverage of the filter by the continuum longward of Lyman-a. 

We have followed the approach of Steidel et al. in calculating 
the effect of this luminosity bias on our sample of z ^ 6 LBGs. We 
compute an effective survey volume using 



Kff(m) = 



dzplni, z) 

az 



where p{m, z) is the probability of detecting a galaxy at redshift z 
and apparent z' magnitude m, and dz ^ is the comoving volume 
per unit solid angle in a slice dz at redshift z. We integrate over 
the magnitude range we are sensitive to, and over the redshift range 
5.6 < z < 7.0 from our colour selection. We calculate that for a 
spectral slope of l3 — —2.0 (i.e., flat in /„) the effective comoving 
volume is 40 per cent the total volume in the range 5.6 < z < 7.0 
(i.e., equivalent to 5.6 < z < 6.1). For our 146arcmin^ survey 
area (after excluding the unreliable regions close to the CCD chip 
gaps) this is a comoving volume of 1.8 x 10^ h^ffMpc"^. For a 
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ID RA Declination z^g i^g («' - ^Oab (^ab ~ -^Vcga) (-'-^fs)vcga Rh SFR^^^.s 

(J2000) (J2000) arcsec Mqjv- 



1 


3'^32™40?00 


-27°48'15'.'0 


25.41 ± 0.07 


26.91 ± 0.22 


1.50 ±0.23 


0.90 ±0.15 


0.72 ± 0.25 


O'.'IO 


17.1 


2 


3h32™34?64 


-27°47'20'.'9 


25.58 ± 0.08 


27.09 ± 0.26 


1.51 ±0.27 


1.94 ±0.10 


1.82 ±0.07 


0'.'18 


14.7 


3 


3h32'"25!59 


-27°55'48'.'4 


24.67 ± 0.03 


26.27 ±0.13 


1.60 ±0.13 






0'.'07 


33.8 


4 


3h32™18!18 


-27°47'46'.'5 


23.71 ±0.01 


25.26 ± 0.05 


1.54 ±0.05 


1.72 ± 0.02 


1.33 ±0.03 


(y.'06 


82.0 


5 


3'^32™38!80 


-27°49'53'.'6 


25.45 ± 0.08 


> 27.3 (3(t) 


> 1.85 {3(t) 


3.12 ± 0.09 


0.22 ± 0.02 


0'.'07 


16.5 


6 


3h32™06?48 


-27°48'46'.'6 


25.35 ± 0.07 


27.17 ±0.28 


1.82 ± 0.29 






0'.'19 


18.1 


7 


3h32™33?19 


_27°39'49'.'2 


25.38 ± 0.07 


> 27.3 (3(t) 


> 1.92 {3 a) 






o'.'n 


17.7 


8 


3h32™23?77 


-27°40'37'.'8 


25.54 ± 0.08 


> 27.3 (3(t) 


> 1.76 (3cr) 






(y.'i5 


15.2 


9 


31^32™ 18! 18 


-27°46'16'.'l 


25.21 ± 0.06 


26.75 ±0.19 


1.54 ±0.20 


1.88 ± 0.11 


2.15 ±0.10 


0'.'19 


20.6 



Table 1. Candidate z > 5.6 galaxies, satisfying our criteria z'^^ < 25.6 mag and {«' — z')ab > 1.5 mag. All magnitudes are apeiture magnitudes within 
a 1 arcsec diameter aperture, and the fluxes have been averaged over the 3 available epochs of ACS data. For those objects that were undetected in i', the 3 a 
limiting magnitude of i^g = 27.3 has been used to place a lower limit on the (i' — z') colour. ' — ' indicates that no data is available. The half-light radius, 
Rh, is defined as the radius enclosing half the flux. Stellar objects have R/^ Ri 0'.'06. Nominal star formation rates (SFRs) are calculated assuming objects are 
placed in the middle of our effective volume (a luminosity- weighted redshift of 2 = 5.8), and are based on the fc-coirected 1500 A rest-frame UV continuum, 
assuming a spectral slope of /3 = —2.0; for a more reddened slope of /3 = —1.1, the star formation rates are 10 per cent higher We use AB magnitudes 
for the HST/ACS images and Vega magnitudes for the near-infrared J and Ks data; the conversions to AB magnitudes are (Jab — ^Voga) = 0-981 and 
(Kab - ^Vega) = 1-843 (Koornneef 1983). 



redder spectral slope of [3 = —1.1 (the average of the z — 3 LBGs, 
Meurer et al- 1997) the effective volume is 36 per cent of the total. 

Since our aim is to investigate the high-redshift progenitors 
of present-day normal galaxies we must survey a comoving vol- 
ume at z ~ 6 sufficient to encompass many ~ L* galaxies by 
the current epoch, in order to reduce the effects of cosmic vari- 
ance. From local surveys, the number-density of L* galaxies is 
0* = (0-0048±0-0006) h^o Mpc"^ (Loveday et al. 1992)- Hence, 
the comoving volume occupied by the progenitors of a present-day 
L* galaxy is V* ~ 200 hjg Mpc''. Our {i' — z') colour selection 
surveys an effective comoving volume of ^ 180, 000 Mpc^. 
Thus in principle, assuming the data reaches sufficient depth, we 
could detect the progenitors of ~ 1000 present-day L* galaxies in 
the CDFS survey field- 
Since significantly fewer objects than this have been identified 
in this work we are able to conclude that at 2 ~ 6 the ancestors of 
L* galaxies have continuous unobscured star formation rates lower 
than ~ 20 h^Q Mq yr^^ per galaxy over 5-6 < z < 6.1. 



4.2 The Star Formation Rate at 2 ~ 6: Lower limits and 
Caveats 

Our selection criteria have identified 6 candidate star-forming 
galaxies with z >5-6, (excluding the probable star, the ERO 
and the potential AGN with nominal SFRs of 17, 15 and 82 
h^Q Mq yr~^ respectively)- These have an integrated star forma- 
tion rate of (122 ± 50) h^Q Mq yr^^ (where the error is based 
on Poisson statistics), giving a comoving volume-averaged star 
formation density between z = 5-6 and z — 6-1 of psfr = 
(6-7 ± 2-7) X 10~* hro Mq yr~^ Mpc~^ 

This value is compared to global star-formation rates found at 
lower redshift in Figure|8| reproduced from Steidel et al- (1999) and 
recalculated for our A cosmology and higher limiting luminosity- It 
appears that the comoving star formation rate (based on the bright 
end of the UV luminosity function) was 6 times less at 2 ~ 6 than 
at 2 = 3 — 4- The interpretation of this result, however, depends on 
several considerations- 

Firstly, the star formation densities plotted have not been cor- 
rected for dust extinction within the galaxies in question. Although 



this will reduce the observed UV flux, and hence the implied star 
formation rate, the method used to select galaxies (the Lyman 
break technique) has been consistently applied at all redshifts above 
2 ~ 2- Thus equivalent populations should be observed at all red- 
shifts and the resulting dust correction will alter the normalization 
of all the data points in the same way if there is no evolution in the 
extinction- 

Secondly, previous authors have integrated the star formation 
density over the luminosity function for their objects down to star 
formation rates of 0- 1 L* even if these SFRs have not been ob- 
served- We prefer to plot what we actually observe and not to ex- 
trapolate below the limit of our observations- If the luminosity func- 
tion is extended down to 0-1 L* then the global SFR is increased 
by a factor of ~ 5 for a = —1-6- 

Furthermore, Lanzetta et al- (2002) note that at high redshifts 
the effects of cosmological surface dimming of objects can lead to 
failure to detect dispersed regions of low star formation and thus 
to an underestimation of the total star formation density when cal- 
culations are based purely on compact and intense regions of high 
star formation- 

Since these potential corrections to our data are expected 
to increase the volume-averaged star formation density, the value 
quoted above should be considered a firm lower limit rather than an 
absolute determination- Combining this work with planned deeper 
HST/ACS observations in the i'- and z'-filters, such as the Ultra 
Deep Field, will enable fainter magnitudes to be probed and allow 
the slope of the luminosity function to be determined (e-g-, Yan et 
al- 2003), further constraining the integrated star formation rate and 
associated contribution to the UV background- 



5 CONCLUSIONS 

We have determined the space density of UV-luminous star-burst 
galaxies at 2 ~ 6 using deep HST ACS SDSS-i' (F775W) and 
SDSS-2' (F850LP) and VLT ISAAC J- and Ks- band imaging of 
the Chandra Deep Field South from the GOODS survey- The Ly- 
man Break technique was used to identify 7 candidate high redshift 
objects, of which one may be an AGN- We reach an unobscured 
star formation rate of ~ 15 h^Q Mq yr~^ at 2 '-^ 6, equivalent to 
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L* for the Lyman break population at z ~ 3 — 4 over a comoving 
volume of ~ 1.8 x 10^ h^Q Mpc^ - which should encompass the 
progenitors of about a thousand L* galaxies at the current epoch. It 
appears that the comoving star formation rate (based on the bright 
end of the UV luminosity function above L*) was 6 times less at 
2: ~ 6 than at z = 3 — 4. These results have used to determine a 
lower bound on the integrated volume averaged global star forma- 
tion rate at jz ~ 6. Hence we have been able to extrapolate the star 
formation history of the Universe to earlier epochs than the exist- 
ing Madau-Lilly diagram. These observations provide an important 
lower limit on contribution of galaxies to the UV background over 
a redshift range that some recent observational evidence suggests 
is the epoch of reionization. 

Note Added in Proof 

We have recently obtained spectroscopy on a subset of the objects 
described in this paper using the DEIMOS instrument on the Keck 
II telescope, in collaboration with Prof. Richard Ellis (Caltech) and 
Dr. Patrick McCarthy (OCIW). While analysis of this data is on- 
going, one of the objects (Object 3) shows Ly-Q in emission at a 
redshift of 5.78, and we report this spectrum in Bunker et al. (2003, 
|astro-ph/0302 4 0'T} submitted to MNRAS (Letters). 
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